^ ; Preprint INFN-NA-IV-94/33 DSF-T-94/33 

a^' 



(N 



> 

(N 

o 



Oh! 



X 



STATISTICS OF PARAMETRICALLY EXCITED 
PHOTON- ADDED COHERENT STATE 

V.V.Dodonov, Ya.A.Korennoy, 

V.I.Man'ko*, and Y.A.Moukhin 

Lebedev Physics Institute, 

Leninsky prospect 53, 117924 Moscow, Russia 

* Lebedev Physics Institute and 

Dipartimento di Scienze Fisiche, 

Universita di Napoli ^^FedericoII", 

Mostra d'Oltremare, Pad. 19 — 80125 Napoli^ Italy 

February 7, 2008 



Abstract 

Photon distribution function, means and dispersions are found explicitly for the nonclassical 
state of light which is created from the photon-added coherent state \a,m) due to a time- 
dependence of the frequency of the electromagnetic field oscillator. Generating function 
for factorial momenta is obtained. The Wigner function and Q-function are constructed 
explicitly for the excited photon-added coherent state of light. Influence of added photons 
on known oscillations of photon distribution function for squeezed light is demonstrated. 



Introduction 



There are several types of nonclassical states of one-mode light, statistics of which has 
been studied last decade. They include squeezed light |]T[], |0, |^, correlated light |^, even 
and odd coherent states P, [^ or Schrodinger cat states 0. All these nonclassical states 
demonstrate either sub-poissonian or super-poissonian statistics of photons. The statistics 
of the coherent photons M , I 



10 



is described by the Poisson distribution function. The 
coherent states of light are considered as the classical states. 

The squeezed light has the property of the reducing the dispersion in one of the pho- 
ton quadrature components. The photon distribution function of squeezed light has wavy 
behaviour for large squeezing [|TT], |]T2|. The quadrature components of the squeezed light 
satisfy the condition of the minimization of the Heisenberg uncertaintly relation |jT3| . 

The correlated states of light 0] satisfy the condition of the minimization of Schrodinger 
uncertaintly relation 
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The photon distribution function of the correlated light has 
the oscillatory behaviour |T^ , |l^ depending on the correlation coefficient of the quadrature 
components which is extra physical parameter characterizing the state of the electromagnetic 
field oscillator. The property of oscillations of the photon distribution function is presented 
also for two-mode light ||18| , [1^ , [^ . 

In Ref.||21|| the state \a,m), which is obtained by the multiple action of the creation 
operator on the usual coherent state has been studied. The photon-added coherent state 
|a,?n) depends on extra discrete parameter, number of excited photons, which infiuences 
the statistics of photons. The statistics of this light including distribution function, means 
and dispersions of quadrature components, Wigner function and Q-function were found 
explicitly. This state of light is interesting because it may be created due to interaction of 

The photon-coherent states of light \a,m) 



one-mode light with the two-level system [21 



have the specific property since they do not contain the contribution of states with fixed 
photon number \n) if the number n is less than m. Thus the distribution function for 
photon-added states p{n) is equal to zero up to the number n = m. 

The photon-added coherent states should be distinguished from the states that are ob- 
tained by the action of the displacement operator on the number state. These states have 



been studied in 
nonzero for n < 



[P^]. For these states the photon distribution function may be 
m. One of mechanisms producing the squeezing and correlation of the 
photon quadrature components is the parametric excitation of the electromagnetic field 
oscillator |2^. The parametric excitation may change the statistics of the photon-added 
coherent states. 

The aim of this work is to study the photon statistics of the photon added states subjected 
to the parametric excitation. We obtain the Wigner and Q-function of the parametrically 
excited photon added states and study the deformations of photon distribution function 
of squeezed and correlated light under the infiuence of the extra physical parameter — the 
quantity of added photons. We will show that the characteristic property of the distribution 
for photon-added coherent states to be equal to zero up to a certain number of photons is 
destroied by the parametric excitation of the electromagnetic field oscillator. 



2 State |a,m) of Stationary Oscillator 

Let us consider the state defined by the formula 



\a,m)=Nexp{ ^)^^ ( e^P(— )l") ) ' (1) 

where N = (m!Lm(— |ap))~^''^ is a normahzation constant, \a) is a coherent state p| and m 
is a positive integer. The state |a) can be defined as 



\a\ 



\a) = exp{ —)exp{aa~^)\0), 



(2) 



where |0) is a vacuum state and a, a"*" are photon creation and annihilation operators. Using 
this series we obtain 

!"•'"> =(m!M-|aP))"^- <" 

This state coincides with the state \a,m) that has been studied in |]21[]. With the aid of (|1]) 
and using the formula for the Hermite polynomial 






de 



(4) 



we can easily write the evident expression for wave function of the state \a, m) in the coor- 
dinate representation 

/I \ 2-"^ H^{q-a/V2){q\a) 

^''"' "^^ = HM-|aP))V2 • (5) 

The states \a,m) should be distinguished from the states \n,a) that are defined as 

\n, a) = D{a)\n). 

These states are related to the shifted harmonic oscillators. The properties and statistics 
of these states were studied in [^, ^, ^]. It easy to notice that the state \n,a) is the 
superposition of the states |a,m) 



\n, a) 



1 



The reciprocal expression can be derived as well 



^ Cr(-«*)"-'"(m!L^(-|ap))V>,m). 



\a, m) = (m!L^(-|an)-i/2 ^ c;;(«*)™-"v^|n, a). 



n=0 



The photon distribution function for the state \a,m) has been found in pT 

\2\ 



p{n) 



n 



!|ap("-"^)exp(-|a| 



[{n — my]'^ Lm{ — \ci\'^)m\ 



(6) 



We could calculate the generating function for this distribution in the form 
which is associated with the expression (^) by the formula 

OD 

G{z) = Y.p{n)z\ (8) 

?i=0 

Below this expression will be generalized for the nonstationary case. 

3 Time— Dependent Photon— Added Coherent States 

We introduce the state |a,m, t) which is ananlogous to the state \a,m) of stationary oscil- 
lator, but the new state satisfies the Schrodinger equation with the Hamiltonian 

2 ^ ' 2 

The initial value for the frequency Q{t) is taken to be ^2(0) = 1. We choose the initial state 
of the parametric oscillator to be photon-added state \a,m, 0) = \a,m), and we will study 
its evolution. In terms of evolution operator U{t) of the oscillator it means that 

\a,m,t) = U{t)\a,m,0)= U{t)\a,m), U{0) = l. (9) 

The evolution operator U{t) is unitary one, 

u{t) u+{t) = i. (10) 

Thus we can write 

\a,m,t) ~ f/(t)a+™ U+{t) U{t)\a) = i+™(t)|a,t), (11) 

where \a,t) is a time-dependent coherent state ||2B|, and 

i+(t) = U{t) a+ f/+(t) (12) 



is the integral of motion of the parametric oscillator. From (^ and (IC) we can derive that 
A'^{t) satisfies the initial condition 

i+(0) = a+. (13) 

By using hermitian conjugation of ([l^) and (|13D we get another integral of motion 

A{t) = U{t) a f/+(t), i(0) = a. (14) 

Due to unitarity property of the evolution operator for hermitian Hamiltonian, the operators 
A{t) and A~^{t) satisfy to the same commutation relation as a and a"*" 

[i(t), i+(t)]= i. (15) 



Thus the algebra defined by the operators A{t), A'^{t), 1 is the same Heisenberg-Weyl 
algebra as in the case of operators a, a+, 1. Using the definition of invariants we get p6| 



-)' ^- R~vM^''"''^ 



Ait) = -^( sit)p - mq), A^t) = --^( e*it)p - e*it)q) (16) 

where c-number function e{t) satisfies the equation 

e{t) + n\t) e{t) = 0, (17) 

with the initial conditions £(0) = 1, e(0) = i, which means that the Wronskian is 

££* - e*e = -2i. (18) 

Thus, the solutions to Scrodinger equation for the arbitrary time-dependence of the fre- 
quency are expressed in terms of the solution e(t) of the classical equation of motion (|1^. 
Let us write down the evident expression for the state \a,m,t) in the coordinate represen- 
tation 

{q\a,m,t) = f r t f \ I2AM/2 ' (^9) 

where {q\(y,t) is the time-dependent coherent state 

iM 1/9 A^q'^ V^aq a^e* laP, ,„„, 

(g|a, t) = T,~^l^e-^l\xv{^ + ^ " ^T - ^)- ^^0) 

We see that the wave function of parametrically excited photon-added coherent state is 
expressed in terms of Hermite polynomials. 

4 First and Second Quadrature Component Moments 

It is easy to derive from (|16D the expressions for photon quadrature components in terms of 
the integral of motion 

p = i=( A{t)e\t) + A\t)e{t)\ q = i=( A{t)e*{t) + A^t)e{t)). (21) 

So, the problem of finding the moments of the quadratures is reduced to finding the average 



values of products of the operators A and ^4+ in the states \a, m, t). Using the results [^ 
we get 



Wi-la 



2) 



(q) = {a, m, t\q\a, m, t) = "'\\' ' {e*a + ea*), (22) 

V2 Lrn{-\a\^) 

{e*a + ea*), (23) 



W(-\a 



V^ Ljni- 



a\ 



(£*V + ^2^*2) L(^)(-|«|2) + 2\e\\m + 1) L„+i(-|«|2) - jep L, 



\a 



2\ 



«^>- " 2L,(-m " " ' " ■ P^' 



.2. _ (e*V + e'a*') L(^\-\a\') + 2\e\\m + 1) L^+i(-|«p) - \e\' L^(-|ap) ,_^ 

^^^" 2L,„(-|an • ^^^^ 

With the help of these expressions it is easy to write dispersions of the quadratures q and 

c^m = m - {q)\ 4(t) = iP") - iP)"- (26) 

All the moments and dispersions are the functions of time as well as of the amplitude |a| 
and the phase of a. In the stationary case the results coincide with the results obtained in 
Ref. m. 



5 Mean and Dispersion of Photon Number 

Let us express the creation and annihilation operators of photons in terms of the integrals 
of motion. We get 

a=^{A{e*+te*)+A+{e + te)), a+ = ^{A+ie - le) + A{e* - te*)). (27) 

The average number of photons in the state \a, m, t) is 

(A^) = (a,m, t|a"'"a|a,m, t) = (|e|^ + |£:|^) 



2 , |.|2n("^ + 1) Lm+i{-\a\'^) 



2 Lmi-\a\ 



I,. ,o ..... 1 Wi-W^ 



^(1^1 +1^1 ) ~ 2 4L {-\oiVV (e + e ) + a (e + £ )), (28) 

and the average of number of photons squared is given by the formula 



lia^af) = {a*\e^ + i^ + a\e*' + e*^) 



^m \ 



\a 



2-\ 



16 Lrr,.(-\a\ 
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4 L„(-|a, , 

+ (m + l)(m + 2){\e^ + ^T + ^e? + 1^1^)^^^^^ 

8 Lm{-\a\^) 

-(m + l)(2(|er + |^r)(kP + l^r + 1) + \e' + e2|2^^-+l(-|«l') 



-(3(|£|^ + \e\') + 2){a*\e' + e') + a\e*' + e*')) 



4L„(-|a|2) 



8Lm{-H^) 



+ i((kl' + 1^1' + 2)' + 2\s' + sY)- (29) 

lo 

With the help of these expressions one can write down the expression for the function which 

describes the dispersion of the number of photons ^7_ 



Q(a,m) = ^^^ . (30) 

It's worth noticing also that the average values of products of any number of operators A" and 
A'^"^ (where m and n are the arbitrary integers) could be evaluated from the commutation 
relations. 



6 Time— Dependent Photon Distribution Function 

Next we study the photon distribution of the field in the state \a,m,t). By definition we 
could write for this distribution fuction the expression 

p{n,t) = |(n|a,m, t)|^. 

To calculate the integral {n\a,m,t) we will use the evident expression for the states \n) and 
I a, m, t) in the coordinate representation, in which the Hermite polynomials are taken in the 
following form 

am 

HUO = lim Q^exp{-(3' + 2^- (31) 

Since f5 is an independent parameter to find the projection of the state |a,m, t) on the 
vector \n) means to calculate the Gaussian integral. Using the Hermite polynomials of two 
variables, that are defined by means of the generating function (see, for example ||25|| ) , 



1 n'^n"^ 

exp{--aRa + aRy)= ^ ^Hi^^y.^y^). (32) 



2 n^O ^'"^' 



we get the distribution function in the form 



p{n,t) 



n!m!(|e|2 + |£|2 + 2)1/2 L„(-|a| 



2^ 



where HJ^^^,!]) is a Hermite polynomial of two variables defined by the 2 x 2— matrix R 
with the following matrix elements 

6 ie* - e 
e* e — IE 



Rn - 


e + ie 
e — ie 


2e 





Rl2 — R2I , ■•Ml 

[e — ie)\e\ 



In the stationary case the matrix R reduces to 



^ -^ (1 oy ' 



and if taking into account that |^| 
where 



/imn = mm[m,n), 



we reproduce the result of Ref. [pll] for stationary case. 



Now we write down the generating function, which is defined as 

oo 

G{z) = J2pin,t)z\ (34) 

n=0 

Calculating the expression for the state \a,m,t) (see above) we derived also a relationship 
between the states \a,t) and \a,m,t). Therefore, it is easy to obtain the relation between 
the projections of these states on the state \n) 

{n\a,m,t) = (m! L^{-\a\^)y^^^exp{ ^)^-^ea;p(— )(n|a, t). 

Let us introduce the distribution function for m = 0. This will be 

Po{n,t) = \{n\a,t)\\ (35) 

Then the relation between the functions p{n,t) and po{n,t) is given by the formula 

p{n,t) = {ml LU-\a\')re-^-^'g-^;^g^^e^-^^o{n,t). (36) 

The same relation exists for the generating functions as well. Due to this to find G{z,t) it 
is enough to compute the function 

oo 

G'{z,t) = Y.Poin,t)z\ (37) 

n=0 

The function po is proportional to the |-ff„P, therefore, we need to obtain the relation for 
summing the modules square of the Hermite polynomials. To obtain it we consider the 
product of two generating functions for the Hermite polinomials 

oo oo ™n m 

exp{-x' + 2xOexp{-y' + 2yC) = Y.Y. -^—MOH^iC), (38) 

making a substitution of variables 

X = \/Xx, y = V\x\ 

where A is a real parameter. Then multiplying both sides of (|38|) by the factor exp{—xx*) 
and calculating the Gaussian integral with the account of the relation 



^ f X*"x""exp{-xx*)dxdx* = n\6„ 



we obtain the expression for G^{z,t), where the relation between A and z is given by the 
following formula 

z\e + ie\ 

A = — j j- 

2\e — ie\ 



Using (|36| ) we find the expression for the generating function 



On 



^(^'^) = -Tr^T\T^.H'^l{o^,o^*)exp[{za, - l)\a\^ - Re{J,,a% (39) 

m\Lm[-\a\'^) 

where J is 2 x 2-niatrix with the following matrix elements 



J 11 — J 22 



ctq 



2e 

J\2 = J2I = —ZCTq, 

4 



One can check that G{l,t) = 1. From the definition (|3^ it follows that 

+ dG{z,t) 

la^a) = hm . 

z~*i dz 

Taking the limit z ^ 1 we see that this result coincides with the expression ( P5| ) for {a^a) 
obtained above. 

7 Quasiprobability Distributions for Field in |a,m,t)— 
State 

Next we study the Q-function which is the diagonal matrix elements of the density operator 
in coherent states basis (see, for example, [0). Thus, the Q-function is defined by formula 

Q{z) = {z\a,m,t){a,m,t\z). (40) 

Let us calculate the matrix element {z\a,m,t). As in the case of finding the distribution 
function the Hermite polynomials can be expressed through the generating function. The 

calculation gives 

exp(-|zp-|a|2)2|^|2-|/7„(^)|2 



Q{z) 



L„(-|a|2)m!(|e|2+|£|2 + 2)V2 
{e + ie)z*'^ + 4:Z*a - [e* - ie 



*^«2 



where 



X exp Re ^^ — , (41) 

\ e — is 



is* \'^'^ „ / 2z* 



i9= , ? = ^ a 

\2ie-ie) ' \e*-ie* 



A transition to the stationary case is obvious. Preexponential factor in the formula (^) 
is related to m added photons and it changes the Gaussian form of the Q-function of 
the coherent state described by the exponent. Thus, tjhe Q-function of the photon-added 
coherent state has the shape of deformed Gaussian with extra maxima and minima. 



The Wigner function of an arbitrarily state, determined by the density matrix p(g, g'), is 
defined by the relation 

+00 

/lY' ly 
p{q H — , g )exp{—ipr)dr. (42) 

—00 
For the pure state \a,m,t) we have 

p{<l,(l') = {q\a,m,t){a,m,t\q'). 



Therefore, to calculate the integral ( ^2]) we make the same computation as for the inte- 
gral {n\a,m,t) found above. For two-dimensional Hermite polynomial having zero diagonal 
matrix elements of the matrix R the following relation holds 

Due to this we can represent the Wigner function in the form 

W{q,p) = I ') ," ';' ' exp{-\tep - lex - v^ap), (43) 

Lm{-\0lY) 

where 7 is defined by 

. . r- /-ee* as* 
7 = ip\e\v2 — iv2——q — — -. 



The physical parameters m describing extra photons produces the factor Lm(|7p) in the 
formula (^31) and this factor changes the Gaussian shape of the Wigner function of the 
coherent state. Since Laguerre polynomial has the wavy behaviour the Wigner function 
of photon-added state has local maxima and minima. Thus, we obtained the phase space 
distributions for the photon-added states of the parametric oscillator. 

8 Numerical Evaluation of Means and Dispersions for 
Specific Time-Dependence of Frequency 

In this Section we consider a spesific time-dependence of the oscillator frequency Vt^{t) in 
form of "the saw with one tooth", i.e., 

fi2(t) = 1, if t0O or t>T, 

and 

fi2(t) = l + kt, if 0<t<T. 

The time- dependence is determined by two parameters T and k. Function e{t) for such a 
"saw" is given by the formula 

£{t) = e'\ if t < 0, 



e{t) = v/rTH(Ci Ji/3(^^^^^) + C^2l^i/3(^^^^^)), ^/ < t < T, 



where Ci and C2 are defined by 



Ci 



■^1/3(3^)^4/3(3^) - J4/3i^)Yi/3i^ 



2 N ' 



^ _ '^i/3(i)(A;-0-'/4/3(i 



|)(fe-0-^4/3(i) 
2 T / 2 AAA / 2 \ 7 . ^ 2 \iA / 2 \ 5 



■^1/3(3^)^4/3(3^) - -''4/3(3fe)>"l/3(3fe 

and if t > T, tlien e(t)-function is given by 



e = Die** + Dse" 



wliere Di and D2 are given by 



g-iT 



/^2 = ^(5(T)+^e(T)). 

At tlie fixed perturbation parameters k and T the dispersions of quadratures q and p will be 
periodical functions of time when t > T. The minimal value of dispersion squared multiplied 
by the factor 2 is the squeezing coefficient s 

1 



H 



We studied numerically the dependence of squeezing on k, a, m at the fixed T = 1. In Fig. 
1 there are presented the results of the numerical calculation for parameter a = 10, which 
corresponds to the multiphoton field. Diferent curves correspond to diferent m. We can 
see that they are very close. Their difference can be seen better at /c ~ 30. In this range 
the values oi m = 1,5, 10,50 correspond to the curves, if one counts from top to bottom. 
The dependence of squeezing on k is an oscillating function. When A; ^ 60 we observe the 
local minimum of squeezing, which is equal approximately to 0.02. In the following local 
minimum the squeezing is even stronger. 

The parameter a = 1 corresponds to the Fig. 2. In this case the field state is intermediate 
between classical and quantum one. We can see that the dependence on m is stronger. In 
the range of k from 30 to 80, the upper curve correspons to m = 5, and the lower one — to 

171 = 1. 

In Fig. 3 there are presented the plots for quantum state of the field a = 0.1. The 
dependence on m becomes even stronger. The curves do not overlap as it was in Figs. 1 and 
2. The number m = 1 corresponds to the upper curve and m = 5 to the lower one. 

So we can conclude that with the decrease of |q;| the field becomes more sensitive to 
photon pumping. 

Now we discuss the behaviour of average number of photons and the dispersions of a 
number of photons. When t > T the Hamiltonian becomes stationary and therefore the 
operator of a number of photons a^a is the integral of motion. The average values of 
(A^), (A^^) and therefore the dispersion of a number of photons do not depend on time when 



t > T. The dependences of the ratio of the dispersion of a number of photons squared and 
mean photon number (A^) (which is denoted Q{a,m)) on k at the fixed T = 1, t > T, a = 1, 
m = 5 is presented in Fig. 4. As in the case of squeezing they have the tendency of the 
increasing role of m when |a| decreases. If we compare Fig. 4 with Fig. 2 we can notice that 
when squeezing increases an average number of photons and the dispersion of it in the field 
also increases. 

Let us discuss the results of the distribution function investigation. For the stationary 
case t > T the distribution p{n) does not depend on time and the most interesting is its 
dependence on n and k. In Fig. 5 p{n) is ploted for m = 2, a = 1. We chose the range of 
k from 50 to 80. According to Fig. 2 it is the range in which strong squeezing is achieved. 
We can see here the known oscillations oi p{n) |Tl|, which proves nonclassical nature of the 
field. At smaller sqeezing these oscillations are smaller, too. The dependence of the photon 
distribution function on k behaves as smoothly oscillating (Fig. 5). While varying k, the 
behaviour of maxima and minima of p{n) is changing. This is seen in Fig. 5 when n varies 
from 15 to 20. 

In Fig. 6 the dependence of the distribution function on a is illustrated. Here we chose 
k = 63, T = 1, 7n = 2, t > T and a is real. At a = we get the distribution function of a 
squeezed Fock state \m,t) [|^ which looks in the coordinate representation as 

m/2 



(g|m,t) = ( 7p ) {m\e\^y^''^exp 




While a is changing in the quantum range from to 1, the distribution oscillations can be 
easily noticed. The increase of a results in the classical limit of a large number of photons. 
It produced a strong delocalization of the distribution function of the parameter n. Thus, 
the oscillations become smaller and invisible in the graph. 

The distribution function atA; = 63,T=l,a = l,t>T for m = 1 and m = 2 is shown 
in Fig. 7. For ?7i = 1 at n > 40 oscillations disappear while for m = 2 they exist at large 
values of n. 

Generaly speaking, the problem of dependence of the distribution function for the squeezed 
and, simalteneously, correlated state on the parameter m requires special investigation. We 
can just suppose that the trend of the dependence on m will be different at different extent 
of squeezing and different correlation between quadrature components. 

In Fig. 8 there is Q-function at T = 1, fc = 10, a = 1, m = 5, t = 10. If a is increased, 
the form of a Q-function peak will resemble the form shown at fig. 8. However, a shape 
and a position will be different. Approximately, the same trend is observed while increasing 
m. At the values of a corresponding to the quantum case the form of peak becomes much 
more sensitive to changes of parameters. It is illustrated in Fig. 9, which represents the 
Q-function at a = 0.1, T = 1, A; = 1, m = 10, t = 10. 

The Wigner function is more sensitive to changes of m. In Fig. 10a there is Wigner 
function at T = 1, /c = 10, a = 2, m = 9, t = 10. The Wigner function of the coherent 
state is the Gaussian function with one maximum. On the other hand, fot the photon-added 
states such shape is deformed and several local maxima are observable. Fig. 10b shows the 
same plot that Fig. 10a but from the direction of p axis. In this figure it is clearly seen that 
the function W{q,p) can be negative. 



In Figs. 11 and 12 the Wigner function is ploted at the same values of parameters 
T = t = l, a = 2,m = l but at different values of k, 10, and 60, respectively. At /c = 10 
the dispersion of q is greater than at A; = 60 but the dispersion of p is less. That is why the 
peak of the Wigner function in Fig. 11 is significantly narrower along q and wider along p, 
than in Fig. 12. 

Thus we demonstrated that parametrical excitation of the electromagnetic field oscillator 
changes the statistics of photon-added coherent states since the squeezing and correlation 
of quadratures produced due to nonstationarity of the oscillator influences the dispersions 
and means of photon numbers in these states. 
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FIGURE CAPTIONS 

Fig. 1. Squeezing of g as a function of k for T = 1, a = 10,, and m = 1, 5, 10, 50. 

Fig. 2. Squeezing of g as a function of k for T = 1, a = 1,, and m = 1, 5. 

Fig. 3. Squeezing of g as a function of k for T = 1, a = 0.1, and m = 1, 5. 

Fig. 4. Tlie dispersion of tlie number of pliotons squared and divided by the average 
number of photons Q{a, m) for T = 1, t > T, a = 1, and 171 = 1, 5, 10. 

Fig. 5. Photon number distribution p{n,t) as a function of n and k for T = 1, t > T, 
a = 1, and m = 2. 

Fig. 6. Photon number distribution p{n,t) as a function of n and a for T = 1, t > T, 
k = 63, and m = 2. 

Fig. 7. Photon number distribution p{n,t) for T = 1, t > T, k = 63, a = 1, and 
m = 1, 2. 

Fig. 8. Q-function Q{z) for T = 1, t = 10, k = 10, a = 1, and m = 5. 

Fig. 9. Q-function g(z) for T = 1, t = 10, A; = 10, a = 0.1, and m = 10. 

Fig. 10a (b). Wigner function W{q,p) for T = 1, t = 10, k = 10, a = 2, and m = 9. 

Fig. 11. Wigner function W{q,p) for T = 1, t = 1, /c = 10, a = 2, and m = 1. 

Fig. 12. Wigner function W^(g,p) for T = 1, t = 1, A; = 60, a = 2, and 771 = 1. 
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